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LARGE DEVIATIONS OE THE EMPIRICAL CURRENT 


IN INTERACTING PARTICLE SYSTEMS 

L. BERTINI, A. DE SOLE, D. GABRIELLI, G. JONA-LASINIO, C. LANDIM 


Abstract. We study current fluctuations in lattice gases in the hydrodynamic 
scaling limit. More precisely, we prove a large deviation principle for the em¬ 
pirical current in the symmetric simple exclusion process with rate functional 
I. We then estimate the asymptotic probability of a fluctuation of the aver¬ 
age current over a large time interval and show that the corresponding rate 
function can be obtained by solving a variational problem for the functional 7. 
For the symmetric simple exclusion process the minimizer is time independent 
so that this variational problem can be reduced to a time independent one. 
On the other hand, for other models the minimizer is time dependent. This 
phenomenon is naturally interpreted as a dynamical phase transition. 


1. Introduction 

In the last twenty years, interacting particle systems have become a main subject 
of research in physics for the insight they provide on the dynamical aspects of sta¬ 
tistical physics. On the mathematical side they provide a source of new interesting 
problems in probability theory. Particularly relevant are the results obtained on 
their hydrodynamical limits and the associated large deviations since a microscopic 
derivation of phenomelogical macroscopic laws can be rigorously established. More 
precisely, for symmetric conservative interacting particle systems, it has been shown 
that the empirical density satisfies a parabolic evolution equation. The associated 
dynamical large deviations rate function measures the asymptotic probability, as 
the number of particles diverges, of fluctuations from the hydrodynamical evolu¬ 
tion. As discussed in , this rate function provides a new approach to the analysis 
of stationary non equilibrium states. These states describe a physical situation in 
which there is a macroscopic flow through the system and the Gibbsian descrip¬ 
tion is not applicable. Rigorous proofs of the dynamical large deviation principle 
have been obtained for some equilibrium models, see e.g. pm, and for the non 
equilibrium simple exclusion process p. 

Beside the empirical density, a very important observable is the current, which 
measures the flux of particle. This quantity gives informations that cannot be re¬ 
covered from the density because from a density trajectory we can determine the 
current trajectory only up to a divergence free vector field. In PP we have intro¬ 
duced, at a heuristic level, the large deviation principle for the empirical current. 
In the present paper we prove it in the case of the symmetric simple exclusion and 
illustrate some relevant applications. 
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The simple exclusion process is a lattice gas with an exclusion principle: a particle 
can move to a neighboring site only if this is empty. The particle dynamics is given 
by a Markov process on the state space { 0 , 1 }^^^^ where = (Z/iVZ)'^ is the 
discrete d-dimensional torus with N‘^ points. We denote by ry = {r]{x),x G a 
configuration, so that r]{x) = 1 when the site x is occupied, r]{x) =0 otherwise. 
Let TT^ the empirical density of particles. The hydrodynamic limit for this model 
is particularly simple: in the limit iV —> oo the empirical density satisfies the 
heat equation. To discuss the large deviation asymptotics we need to introduce the 
mobility x(7r), which describes the response to an external held; for the symmetric 
simple exclusion we have x(7r) = 7r(l — tt). We introduce the integrated empirical 
current W^, which measures the total net how of particles in the time interval 
[0, t], associated to a trajectory t ]. We shall prove a large deviation principle which 
can be informally written as follows. Fix a possible path Wt, t G [0,T] of the 
integrated empirical current, then 

P^(Wf « Wt , t G [0,T]) ~ exp { - iV''/[o.T](W)} (1.1) 

where the rate functional is 

+ + ( 1 . 2 ) 

In the above formula Wt is the instantaneous current at time t. Moreover TTt, which 
represents the associated huctuation of the empirical density, is obtained from W 
by solving the continuity equation i9t7r-|-V-W = 0. Finally (•, •) denotes integration 
with respect to the space variables. Note that equation can be interpreted, 
in analogy to the classical Ohm’s law, as the total energy dissipated in the time 
interval [0,r] by the extra current W + |V7r. The large deviation principle of the 
empirical density m, as we show, can also be easily deduced from (n~T]i-in~a. 

Using fl~Tll-01l we then analyze the fluctuations properties of the mean em¬ 
pirical current W^/T = T~^ dt'Wt over a large time interval [0,T]. This is 
the question addressed in |H] in one space dimension by postulating an “additivity 
principle” which relates the fluctuation of the current in the whole system to the 
fluctuations in subsystems. We show that the probability of observing a given, 
divergence free, time averaged fluctuation J can be described by a rate functional 
$(J), i.e. as iV —» oo and T ^ oo we have 

pAT ^ j) _ gxp { - T<i>(J)} 

The functional $ is characterized by a variational problem for the functional d[o,T] 

$(J)= hm infl/[o.T](W) (1.3) 

1 —W 1 

where the infimum is carried over all paths W^ such that Wy = TJ. 

Let us denote by U the functional obtained by restricting the infimum in 
to the paths W such that Wt is divergence free for any t G [0,T]. The associated 
density profile TTt does not evolve. We get 

um =i„f 1([J + + Iv,]) (1.4) 

where the infimum is carried out over all the density prohles p. This is the functional 
introduced in in one space dimension. For the symmetric simple exclusion 
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process we prove that the additivity principle postulated in [H] gives the correct 
answer, that is $ = ?7. On the other hand, while <i> is always convex the functional 
U may be non convex. In general Oil we interpret the strict inequality $ < ?7, 
as a dynamical phase transition. In such a case the minimizers for OS become 
in fact time dependent and the invariance under time shifts is broken. In ^ we 
have shown that, for the one-dimensional Kipnis-Marchioro-Presutti (KMP) model 
0 cn, which is defined by a harmonic chain with random exchange of energy 
between neighboring oscillators, the following holds when it is considered with 
periodic boundary condition. The functional U is given by 17( J) = (1/2) J^/x(m) = 
{lf2)j‘^f'm?, where m is the (conserved) total energy. Moreover, for J large enough, 
$( J) < U{J). This inequality is obtained by constructing a suitable traveling wave 
current path whose cost is less than U{J). In the present paper we give a more 
formal presentation of these results. We finally mention that the strict inequality 
$ < [/ also occurs for the weakly asymmetric exclusion process for sufficiently large 
external held [7|. 

2. Notation and results 

For N > 1, let = (Z/NZ)"^ be the discrete d-dimensional torus with N‘^ 
points. Consider the symmetric simple exclusion process on T^. This is the Markov 
process on the state space Xn ■= {0,1}^" whose generator Ljv is given by 




\x-y\ = l 


In this formula, rj = {rj{x), x G T^} £ Xjy is a conhguration of particles, so that 
ry(x) = 0, resp. ry(x) = 1, if and only if site x is vacant, resp. occupied, for rj, and 
is the conhguration obtained from rj by exchanging the occupation variables 
rj{x), r]{y): 



{a-^'yTj){z) = < T]{y) ifz = x, 
[ rj{x) li z = y . 


Notice that we speeded up time by . Denote by {??( : t > 0} the Markov process 
with generator Lm and by P|))v the probability measure on the Skorohod space 
, Xjy) induced by the Markov process ijt and a probability measure jj,^ on 
Xfr, standing for the initial distribution. When is a Dirac measure concentrated 
on a conhguration , we denote by . An elementary computation shows 
that this process is reversible with respect to any Bernoulli product measure on 
{0, with parameter m £ [0,1]. 

Denote by the space of hnite signed measures on the d dimensional 

torus of side 1, endowed with the weak topology and hy T = d>x(T‘^) the set of 
positive, measurable functions bounded by 1 endowed with the same weak topology. 
For a hnite signed measure m we let (m, F) be the integral of a continuous function 
F : » K with respect to m. Likewise for a prohle tt £ F and F £ C'(T‘^) we 

denote by (tt, F) the integral of ttF. 
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For TV > 1 and a configuration 77 S Xn, denote by = Tr^{r]^) the empirical 
density of particles. It is defined as 

TT^(u) = ^ 1{mTV g B{x)}ri{x) , 

where, for x = {xi ,..., Xd) & B{x) is the cube [xi,xi + 1) x • • ■ x [xd, Xd + !)• 
Set TT^ = TT^(rjt) and notice that belongs to J- for each f > 0. 

For t > 0 and two neighboring sites x,y G Tff, denote by the total number 
of particles that jumped from a; to y in the macroscopic time interval [0, t]. Let {e^ : 
1 < A: < d} be the canonical basis of the difference _ jx+ej,x 

represents the net flow of particles across the bond {x, x + Cj} in the time interval 

[0,t]. 

For t > 0, we define the empirical integrated current = {W ^^,..., g 

Md = as the vector-valued finite signed measure on induced by the 

net flow of particles in the time interval [ 0 ,t]: 

WN ^ ^-(d+ 1 ) ^ ^ j = I,..., d 

where Su stands for the Dirac measure concentrated on u. Notice the extra factor 
N~^ in the normalizing constant which corresponds to the diffusive rescaling of 
time. In particular, for a continuous vector held F = {Fi, ..., Fd) g C'(T‘^;M‘^) the 
integral of F with respect to Wf'^, also denoted by (W^, F), is given by 


d 

(Wf,F) = TV-(^+i)^ . (2.1) 


The purpose of this article is to prove a large deviations principle for the empirical 
integrated current and discuss the asymptotic behavior as t —> 00 . We start 
with the law of large numbers. Fix a profile X G F and let : TV > 1} a sequence 
of measures on associated to A in the sense that the empirical density converges 
to A in probability with respect to . Namely, for each F G C(T‘^) and d > 0, we 
have 


lim 

N—^CXi 



Ijd 


X{u)F{u)du >1^1 


= 0 


( 2 . 2 ) 


It is well known, see e.g. cni, that in such a case the empirical density converges 
in probability to p = p{t, u) which solves the heat equation 


/ dtp = (l/2)Ap, 
I P(0,-) = A(-), 


(2.3) 


where A = V• V stands for the Laplacian and V for the gradient. We claim that the 
empirical current converges to the time integral of —(l/2)Vp. This is the content 
of the next result which is proven in Subsection EH in a more general context. 


Proposition 2.1. Fix a profile X G F and consider a sequence of probability mea¬ 
sures associated to X in the sense of (EH). Let p be the solution of the heat 
equation 12.311 . Then, for each T > 0, <5 > 0 and F g C'(T'^; R'^), 


lim P 

N^OO 


N 


(W. 


T,r) 


+ ( 1 / 2 ) dt F{u) ■ Wp{t,u) du 


/■fd 


> S 


= 0 . 
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We turn now to the large deviations of the pair Fix T > 0 and 

denote by D{[0, T],Md x T) the space of cad-lag trajectories with values in Md x T 
endowed with the Skorohod topology. Fix a profile 7 S C^(T‘^) bounded away from 
0 and 1: there exists 5 > 0 such that 5 < 7 < 1 — 5. To focus on the dynamical 
fluctuations, we assume that the process starts from a deterministic initial condition 
which is associated to 7 in the sense that 71^(77^) ^ 7 in iF. 

Let be the set of trajectories (W, tt) in D{[Q,T]^ Md x J-) such that for any 
t e [0,r] and any F e C\T‘^) 


(7rt,F)-(7,F) = (Wt,VF), Wq = 0 (2.4) 


Note that Slj is a dosed and convex subset of (^([O, T], A4d x Equation (I2.4|l is 
the weak formulation of the continuity equation -|- V • Wj = 0, tiq = 7 , Wq = 0. 
For each F € C^’^([0,T] x define the convex and lower semi-continuous 

functional Jp as follows. If (W,7r) £ C{[0,T]-,Md x F) we set 

Jf(W,t) = (Wt,Ft) - [ dt{Wt,dtFt) 

Jo 

- dt{xM,\Ft\^) , (2.5) 

where, here and below, x(a) = a(l — a) is the mobility. We set JF(W,7r) = -too 
if (W,7r) ^ C([0,T]; A4£; x F). Since C{[0,T];Md x F) is a, closed subset of 
D{[0,T]-,Md X F), the argument in ^1 §10.1] proves the convexity and lower 
semi-continuity of Jp. 

Let finally 


J(W,7r)= sup JF(W,7r), /(W, tt) 

Feed 


J(w,7r) if(W,T)Ga^, 
-too otherwise. 


Notice that the functional J is convex and lower semi-continuous, properties 
which are inherited by I because 2l.y is a closed and convex subset of (^([O, T],Jd[^ x 
F). Notice furthermore that the continuity equation (12.411 determines the trajectory 
TT as a function of W and the initial condition 7. Therefore, the rate function 
/(W,7r) can be thought as a function of W and the initial density profile. In 
Subsection El we derive a more explicit formula for the rate function EH). If 
/(W,7r) < 00 we have 

= \ dt{x(^,),\F,\^) , 

where the vector-valued function Ft is the solution of 


dtWt + (l/2)V7rt = x(^t)Ft . 

Thus, formally, 

where Wf = dtWt is the instantaneous current at time t for the path (W,7r) and 
Wt (tt) = —(l/ 2 )V 7rf is the typical instantaneous current at time t associated to 
the density profile Tit. Recall in fact that the hydrodynamic equation (lOl is in our 
case the heat equation. 
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Theorem 2.2. Consider a profile 7 G C'^(T‘^) bounded away from 0 and 1 and a 
sequence {rj^ ■ N > 1} such that {rj^) ^ j in tF. Then, for each closed set F 
and each open set G of D{[0,T], A4d x F), we have 


Mm —\ogr\ {W^\F^)€F <- inf 

N^oa Nd <=> V J (W,7r)eF 


Jim ^iogP71(W",.")eG] > 


We next state a “large deviation principle” for the mean empirical current /T 
in the interval [0, T] as we let first N ^ 00 and then T —>■ 00 . 

Let us denote by ;B C Md the set of divergence free measures, i.e. 

S := {j G Xd : (J, V/) = 0 for any / G (2.7) 

which is a closed subspace of Md- Given m G (0,1), we introduce the set of profile 
with mass m, i.e. we set Fm '■= {p & F : Jjddu p{u) = m}. We finally define 
Urn. ■ Md [ 0 , + 00 ] as 

Um{3):= inf ^([j + ^Vp], ^ + (2-8) 

2 \ 2 p(l - p) 2 / 

0<p<l 

if J G , J(d'u) = jdu, and 17m(J) = +00 otherwise. In Section I^TI we show that 
Um is a lower semi-continuous convex functional. 

In Section o we also prove that in the one dimensional case, where J € B, 
J{du) = jdu, implies that j is constant du-a.e., we simply have 

(1 

, , -;- 7 if J = jdu for some j du-a.e. constant 

U^{J) := <^ 2 m(l - m) (2.9) 

I -|-oo otherwise 


Theorem 2.3. Let m G (0,1), 7 G C'^(T‘^) n Fm bounded away from 0 and 1, and 
T]^ G T/v a sequence such that {rj^) ^ j in F. Then, for each closed set C and 
each open set G of Md, we have 


lim lim 

T—>00 N — >CX. 




G G 


lim lim TfTjd 

T-,00 N^oo d Jy 


I 

- W-T. 

IT ^ 
I 

IT 


< 


logP^« 


-W^ G G 


- inf C/m(J) , 

J GO 

> - inf c/m (J) ■ 


3. Large deviation for the empirical current on a fixed time interval 

In this Section we prove Theorem E3 The proof is similar to the one of the 
large deviations principle for the empirical density, see m or [TT)1 Chapter 10]. We 
therefore present only the main modifications. 

3.1. Weakly asymmetric exclusion processes. In this subsection we prove the 
law of large numbers for the empirical current of weakly asymmetric exclusion 
processes. Proposition o follows as a particular case. 

Fix T > 0 and a time dependent vector-valued function F = {Fi, ..., Fd) G 
C'^’^([0,T] X T'*;IR'^). Denote by Lf,n the time-dependent generator on Xn given 
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by 


iV^ 

{L^.Nf){v) = — 


d 

EE 

i=l a:GT^ 


^x,x-\-ej 




where the rate (t, 17 ) is given by 

ry(a;)[l - r){x + e,)] + r){x + ej)[l - v{x)] . 

Hence, for N large, instead of jumping from x to x + ej (resp. x + Cj to x) with rate 
1/2, at a macroscopic time t particles jump with rate (1/2){1 + N~^Fj{t,x/N)} 
(resp. (1/2){1 — N~^Fj{t^x/N)]) and a small drift appears due to the external 
field F. For a probability measure on Xjq, denote by Pp the measure on the 
path space Z?(R+, Xjsi) induced by the Markov process rjt with generator Ly,n and 
initial distribution /i^. 

Let be the unique weak solution of the parabolic equation 

r = (l/2)Ap-V-{x(p)F}, 

I P(0,-) = A(-). 

Write the previous differential equation as 

pt + V • WF(pt) = 0 , 


where Wf(p) is the instantaneous current associated to the profile p and is given 

by 


Wf(p) = -(l/2)Vp + x(p)F . 

The main result of this section states that 'Wf converges in probability to the time 
integral of Wf(p): 


Lemma 3.1. Fix a profile A : » [0,1] and consider a sequence of probability 

measures : N > 1} on Xjy associated to X in the sense of (12.211 . For each 
t > 0, 5 > 0, G e G(T'^) and H G 


N—^oo 

lim P£(, 

AT^oo 


7rf,G) - (pP,G) 


> (5 


= 0 


(Wf,H) - / ds(WF(pr^),H) 


> 5 


= 0 , 


where (WF(pf’^),H) stand for 


(l/2)(p^^V•H} + (x(p^'^),F,•H). 


Proof. The law of large numbers for the empirical density follows from the usual 
entropy method, see e.g. cni Chapter 6 ]. For each t > 0 the empirical density 
converges in probability to p^’^(t, •)• 

To derive the hydrodynamic equation for the current, fix t > 0 and a smooth 
vector field H : ^ Let be the martingale defined by 

__ rt Ar2 

Wf.H = (Wf,H)- / + 

i,x 

ft at2 

+ I ^^ 2Nd+-t E . 












L. BERTINI, A. DE SOLE, D. GABRIELLI, G. JONA-LASINIO, C. LANDIM 


An elementary computation shows that the quadratic variation of this martingale 
vanishes in L^(Pp ^jv) as t c». On the other hand, after a Taylor expansion and 
few summations by parts, the time integral can be rewritten as 

ds(^f,(l/2)V-H + F,-H) + Of,h(A'-') 

“ / ^ Hji^/^)^As^x/N)r]s{x)T]six + ej) , 

i=i rceT^ 

where Of,h(-^~^) is an expression whose absolute value is bounded by CN~^ for 
some constant depending only on F and H. By the two block estimates and the 
law of large numbers for the empirical density, as TV t oo, the previous expression 
converges in ^^-probability to 

r ,/i. 

{1/2) f ds{p^’\V-U) + rds(x(pP),F,.H) . 

Jo Jo 

Since the martingale vanishes in as TV t oo, the lemma is proven. □ 

The same result holds for the generator Lf,a defined by 

Af2 

{LF,Nf){v) = — Y. 

\x-y\ = l 

However, the computations of the exponential martingales in next section are 
slightly more complicated if we use this expression instead of Tf,a- 



3.2. Large deviations upper bound. We first remark that Lemma ESI below 
implies that the probability of the event (W^,7r^) ^ C'([0,T];7Wd x IF) is super- 
exponentially small as ^ oo. Recalling the definition EHl) of the rate function I, 
it is therefore enough to prove the upper bound for closed subsets of (^([0, T]-,Md x 
J^). We shall first prove it for compacts and then show the exponential tightness. 

We start by recalling the super-exponential estimate of Pin]. For a positive 
integer £ and x in Z'^, denote by rf{x) the empirical density of particles on a box 
of size 2^ -I- 1 centered at x: 


r]\x) = 


{2e- 


l)d ^ 
\v-x\<e 


v{y) ■ 


Moreover, for 1 < j < d, e > 0, let 




1 


E 


1 

(2£ivTTF 


E d{y)viy + e-j) - [v^^x)]^ 

\y—x\<eN 


Theorem 3.2. For each 1 < j < d, T > 0, each sequence of measures : N > 
1} and each d > 0, 


lim sup lim sup log P^w 


Vj,N,e{Vt) dt 


> s 


—oo . 
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Fix a vector-valued function F : [0,T] x ^ in and denote by 
dFp /dP^^iT) the Radon-Nikodym derivative of with respect to P^w re¬ 

stricted to the time interval [0,T]. A long but elementary computation gives that 


1 dP^ 1 ^ fT 

E [ m-m-iwr*-’ 


Nd 


dF\ 


i=i xeT %' 


rp 'J' d 

^ dt(7rf,V-F*)-i^ ^ F,(t,cr/lV)V,/i,(r;0 + OF(iV-i), 


i=l xGTi 


where hj(r]) = r]{0) + rj{ej) — 2rj{0)r]{ej) and denotes the translation of x. In 
particular, on the set 


d 

El 

i=i 


Vj,N,s{vt)dt 


< s, 


integrating by parts the first term on the right hand side of the penultimate formula, 
we obtain that logjdPp 


,/dF^!^}(T) is bounded below by 


- C(F){e + <5} 


for every e > 0. Here, 


JF,e.5(W,T) 


{Wt,Ft) - r dt{Wt,dtFt) 

Jo 


C'(F) is a finite constant depending only on F and 7rf is the function defined by 
TT^{u) = J[u-si u+eX\ dv'Xt{v), where [u — el,u + el] is the hyper-cube [iti — 

e, Ml + e] X • • ■ X [ud — e, Md + ej- 

For each e > 0, d > 0 and F of class C^d, the functional Jf,e,s is continuous in 
C([0,r]; Add X T). By repeating the arguments presented in^| § 10.4 ], we then 
obtain that for each compact set K of C([0,r], Add x T), 


limsup^logP 

N —>-oo dy 


N 


iW^,7r^)GK 


< 


- inf J(W, tt) 
(W.7 i-)gA' 


where J is defined in (ESI). 


To extend the upper bound from compact to closed sets, we next prove the 
exponential tightness of the sequence (W^,7r^). As stated before, the following 
Lemma also implies that the probability of the event (W^,7r^) ^ (^([O, Tj; AddX A") 
is super-exponentially small. 


Lemma 3.3. Fix a sequence of measures : N > 1}, a continuous function 
G : —> R, a vector-valued function H : ^ R"^ in and £ > 0. Then, 


lim limsup logP^^w 

5^0 N^oo Nd ^ A* 


sup 

*- |t —s|<5 


< TT; 


^ G > - < 7rf,G > 


> e 


= —oo 


lim lim sup ^ log Pl^V 

5-.0 Nd ^ 


sup I 


< Wf ,H > - < Wf,H > I > £ 


= —CX) . 
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Proof. The proof of the first estimate is similar to the one of the exponential tight¬ 
ness of the empirical measure presented in nm § 10.4]. In our context the initial 
configuration is not however the invariant measure. The necessary modifications 
are worked out below in the proof of the second statement of the lemma. 

To prove the second estimate, first observe that by a triangular inequality and 
since 


limsupA^ '^logjoAr-I-< max< limsup log oat , lim sup ‘^log&Afl 

iV—>-oo ^ ^^ 


N—*oo 


it is enough to estimate 
max lim sup log 

0<fe<T5-i at^oo A^"^ ^ 


sup 


< wf, 


N- 




(3.2) 


> e/3 


(3.3) 

where tk = kd. By (E3, we may also disregard the absolute value in the previous 
expression provided we estimate the same term with —H in place of H. Fix a > 
0 and denote by Mt = A4t(a,H) the mean one exponential martingale whose 
logarithm is given by 


a 

N 






ds ??s(a:)[l 





ds r^s{x + e,)[l - 



Since Hj are functions, a Taylor expansion and a summation by parts show that 
the expressions inside the integrals in the last two terms are bounded by Chu)! -|- 
a)N‘^, where Ch is a finite constant depending only on H. Therefore, by multiplying 
by aN‘^, adding and subtracting the appropriate integrals and exponentiating, we 
get that 




< . 
— 


sup <Wf,H>-<W",H>> £/3 
^ sup Mt/Mtk > exp{{l/6)aN‘^e} 

tk ^t‘^tk-\-l 


provided Ch)! + a)S < e/6. Since is a positive martingale equal to 1 

at time tk, by Doob’s inequality, last expression is bounded by exp{ —(l/6)aA^'^e}. 
Therefore, (lOl is less than or equal to —ae/6 for all S small enough. This shows 
that the second expression in the statement of the lemma is bounded by —ae/6. 
Letting a / cx), we conclude the proof. □ 


Standard arguments, presented in d § 10.4], together with Lemma [3.31 permit 
to extend the upper bound for compact sets to closed sets. 

We conclude this section proving that we may set J(W,7r) = -|-cx) on the set of 
paths (W,7r) which do not belong to 
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Lemma 3.4. Fix a sequence of measures ■ N > 1} and a function H : 

[0,T] X T"* m Let 

- {n^,Ho) - [ dt {irf, dtHt) , (3.4) 

JQ 

^ f {d^^H){t,x/N)dWr+^^ . 

i=i xGf% ° 


Then, for any 5 > 0, 


l™sup^logP^„ 


VTi7T,W^\H) 


> 5 


—oo . 


Proof. Fix a function FI : [0,T] x T'^ ^ K in A summation by parts shows 

that 


1 

Jfd+I 



N{H{t,x + ej/N) - H{t,x/N)}dWf'''"^''^ 


1 


^ / H{t,x/N)Y,d{W: 

xgtI i=l 


- W, 


x,x-i-ej 


}■ 


Since ~ increases by one each time a particle jumps to x 

and decreases by one each time a particle leaves x, this sum is equal to r]s{x) — r]o{x). 
In particular, an integration by parts, gives that the previous integral is equal to 
,H) defined in 13.411 . Therefore, by a second order Taylor expansion, 


VT{7r^,W^,H) = ^ f {dlM){t,x/N)dW:’- 

i=i xgt^ ° 

, OLr(l) r T®.a;+ej , jx+ej,x^ 

^ ]^d+2 2^ 2^ J ’ 

i=l x^Tf 


where oh{I) depends on H and vanishes as A^ t oo. 

We prove that the first expression on the right hand side is super-exponentially 
small, the argument for the second one being similar. To keep notation simple, let 
Fj = df.H. By Chebyshev inequality, for every a > 0, 


t,N 

■m" 


]\[d+2 


J2J2 F,{t,x/N)dWf’ 


X + Si 


1=1 xeTf,' 


> S 


< 


(3.5) 


— aSIV^joN 

e Ib^jv 


exp{a7V-2|^ ^ F,{t,x/N)dWr^^^\}] . 

•_ 1 _rm^ 0 


1=1 xeTi 
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Since + e ^, we estimate this last expectation without the absolute value. 

Denote by Mt the mean one exponential martingale whose logarithm is given by 

E r F,it,x/N)dwr^^^ 

i=i xeT% ° 

1 = 1 xGTff ° 


-iV^E E r dtT^tix + e,)[l - rjtix)]{e-‘^^~"Fit,-/N) _ I 

1=1 xeT% ° 

Since Fj are continuous functions, a Taylor expansion and a summation by parts 
show that the last two integrals can be written as 

1 = 1 xGT% 


where of(1) is an expression depending on F which vanishes as iV | oo. 

Since is a mean one exponential martingale, the right hand side of (13.511 is 
bounded above by 


expafV'^l-(5 + OF(l)dr + C(F)a(iiV-2T| . 


In particular, 


1 


N —>-oo d.\ 


1 


]\J-d+2 


EE/ dtFj{t,x/N)dWt 


x,x-\-ej 


1=1 xeT% 

for every a > 0. This proves the lemma. 


> 5 


< —aS 


□ 


From this lemma, (1^ . and Lemma lO it follows that for every closed set F, 
every ^ > 0 and every finite family {Hj, 1 < j < £} oi functions in 


hm sup ^ log , TT^) e F 

N —>-CXD 


< - inf J(W, tt) , 

(w.7r)eFnai 


where 


t 

2 lf = n {(».») : |FT( 7 r,W,iLd)| < . 


1=1 

Since this inequality holds for every (5 > 0 and every finite sequence FIj, letting 
<5 I 0 and considering a dense family of functions Hj, we obtain that 


1 


limsup-^logPV (W\7r''')GF <- inf J(W,7r), 

N‘^ ^ ' i - (w.77)6Fna ' ' 


where 21 is the set of paths (W,7r) such that lit + V • Wj = 0. Up to this point, 
we did not need any assumption on the sequence of initial measures /i^; but the 
hypothesis that we are starting from a deterministic profile plays now a role to 
replace the set 21 in the previous formula by the set 21.^, proving the upper bound 
of the large deviations principle. 
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3.3. The rate function. To prove the lower bound of the large deviations principle 
in Theorem 12.21 we first obtain an explicit representation of the functional / on the 
paths with finite rate function. 

Given a path tt £ £>([0, T]; IT), we denote by the Hilbert space of vector¬ 

valued functions G : [0, T] x ^ endowed with the inner product (•, •),r defined 

by 

(H, G} 7 r = / dt f dux{'^{t,u))tl{t,u) ■ G{t,u) . 

Jo JT'i 

Fixapair (W, tt) such that /(W,7r) < oo. In particular (W,7r) G T]; TWd x 
£”). Following the arguments in na § 10.5], from Riesz representation theorem, we 
derive the existence of a function G in such that 

/(W,7r) = d<(x(7rt), lGt|2) 

dtWt + (l/2)V7r, = xMGt . 

where the last equation has to be understood in the weak sense: for each H £ 
K'^) and each 0 < s <t <T, we have 

(Wt,H) - (W«,H) = (1/2) f dr(7r„V-H) -b 


dr(x(7r,),G,-H) 


3.4. The lower bound. In this Subsection we prove the lower bound in Theo¬ 
rem E21 Denote by S the set of trajectories (W, tt) in 21.^ for which there exists 
a vector-valued function G in C'^’^([0,T] x T'^) such that (W,7r) is the solution of 

(1^ . 

For paths (W, tt) in S, we may repeat the arguments presented in the proof of 
the lower bound in the large deviations principle for the empirical density in uni 
§ 10.5] to conclude that for each open set G 


lim inf logP^)^ 


(W 


N 


)£G 


> - inf ICW, tt) . 
(W.7r)eGn5 


To conclude the proof, it remains to show that for all pairs (W, tt) with finite rate 
function, /(W, tt) < oo, there exists a sequence (W^, tt^,) in S converging to (W, tt) 
and such that lim^^oo d(Wfe, tt^) = /(W, tt). When this occurs, we shall say that 
the sequence (Wk,7Tk) /-converges to (W, tt). In the context of the symmetric 
exclusion process, the argument is not too difficult because the rate function is 
convex. We follow the proof of the lower bound presented in [2]. 

The proof is divided in two steps. We first show that there exists a sequence 
(Wfc,7rfc) which /-converges to (W,7r) and such that, for each k, nk is bounded 
away from 0 and 1 uniformly in [0, T] x To do it, following |2], we consider 
a convex combination of (W, tt) with the solution of the hydrodynamic equation 
EU with external field F = 0 and initial condition ttq = 7 , Wq = 0. 

Consider now a pair (W, tt) whose empirical density tt is bounded away from 
0 and 1. Since /(W,7r) is finite, by Subsection 13.31 there exists a vector-valued 
function G in LJ{tt) satisfying (13. till . Since tt is bounded away from 0 and 1, L‘^{tt) 
coincides with the usual LJ space associated to the Lebesgue measure on [0, T] x T'^. 
Consider a sequence of smooth vector-valued functions G„ : [0,T] x ^ 
converging in LJ to G and denote by (W",7r") the pair in 21^ which solves (13.611 
with Gn instead of G. Repeating the arguments presented in 0 § 3.6] , one can 
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prove that (W",7r") /-converges to (W,7r). This concludes the proof of the lower 
bound. 


3.5. Large deviations for the empirical density. In this subsection we show 
that the large deviation principle for the empirical density, proven in na, follows 
from Theorem Indeed, the large deviations principle for the empirical density 
can be recovered from the one for the current by the contraction principle. The 
rate function X is given by the variational formula 

X{tt) = inf /(W,7r), (3.7) 


where SUtt stands for set of currents W satisfying lit -|- V • Wt = 0, as formulated 
in 112.4|l . 

This variational problem is simple to solve. Let us first assume that tt is smooth 
and bounded away from 0 and I. Fix a current W in 211,^ and denote by G the 
external field associated to W through equation (13.till . For 0 < t < T, let Ht be 
the solution of the elliptic equation 

V • (x(7rOGt) = V • {xMVHt) 

and set Ft = x(7rt){Gt — V//t}. By definition V • Ft = 0. Let w be the current 
defined by 

Wt + (l/2)V7rt = x(7rt)Vi/t . (3.8) 

w belongs to 2U,r because by construction V-w = V• W. Moreover, by the explicit 
formula (1^ for the rate function and by definition of F, 

I{W,7t) = dt{xint),\Gt\^) 

= 11^ dt{{xM, |V//tp) + 2(Ft • Vi/t) + (x(7rt)-i|Ftp)} . 

Since V • F = 0, an integration by parts shows that the cross term vanishes. On 
the other hand, by the explicit formula for the rate function and by (13.811 . the 
first term on the right hand side is /(w,7r). Thus, 


/(W, tt) > /(w, tt) . 

In particular, in the variational problem , we can restrict our attention to 
currents W for which the associated external fields G are in gradient form. 

Now, consider two currents W^, in SUtt and assume that both external fields 
Gi, G 2 associated to these currents through 113. till are in gradient form: G^ = V//^ , 
j = 1, 2. Taking the divergence of (1^ and recalling that Tft -f V • Wt = 0, we 
obtain that 

it = (l/2)A^t - V • ixM^m) 

for j = I, 2 and each 0 <t <T. In particular, V • [xi'^t)'^[Ht ~ Ht]) = 0- Taking 
the inner product with respect to — Hf and integrating by parts, we get that 



dux(7rt)|V//t^ 


= 0 


for every 0 < t < T. In particular, /(W^,7r) = /(W^,7r). This proves that the 
variational problem (TTtIi is attained on currents for which the associated external 
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field is in gradient form: 

X( 7 r) = inf IiW,7T) = i r dt{x{7rt)\VHt\^) , 

wean, z Jq 

where Ht is given by 

TTt = il/2)A7Tt - V • {xMVHt) . 

This is exactly the large deviations rate function for the empirical density obtained 
in ini This identity has been obtained for smooth paths tt bounded away from 0 
and 1. However, by the arguments of the previous subsection, we can extend it to 
all paths TT. 


4. Large deviations of the mean current on a long time interval 


In this Section we investigate the large deviations properties of the mean em¬ 
pirical current 'W^/T as we let first —> oo and then T oo. We emphasize 
that the analysis carried out in this section does not depend on the details of the 
symmetric simple exclusion process so that it holds in a general setting. 

Given a profile 7 G C'^(T'^), T > 0, and W G I?([0, T]-,Md), let tt G T*([0, 
be the solution of (12.411 and denote by .f[o,T](W| 7 ) the functional defined in 112.till , 
in which we made explicit the dependence on the time interval [0,r] and on the 
initial profile 7 . We define 4 ) 7 ’(-17) : —> [0, -l-oo] as the functional 

‘fT(J|7) = ^^mf /[o,T](W|7) (4.1) 

1 WgAt.j 


where 

At,j := {W G D{[0,T]-Md) : Wt = Tj} 


Recalling that the set B of divergence free measures has been defined in (E3 we 
also define 


4>(J|7) := 


f inf 
I T>0 

[- 1-00 


if J G H 
otherwise 


(4.2) 


Denote finally by 4 )(J| 7 ) := sup^gj infj/g;/ 4 >(J'| 7 ), where U C Md is open, the 
lower semi-continuous envelope of 4 >(-| 7 ). 


Remark 4.1. The functional 4 >(-| 7 ) depends on the initial condition 7 only through 
its total mass m = f dwjiu)- This holds in the present setting of periodic boundary 
conditions; in the case of Dirichlet boundary conditions, when the density is fixed 
at the boundary, 4 >(-| 7 ) would be completely independent on 7 . Furthermore the 
functional 4 )(-| 7 ) is convex. 


Theorem 4.2. Let 7 G C'^(T‘^) bounded away from 0 and 1, and G Xn a 
sequence such that ir ^^ 7 m .F. Then, for each closed set C and every open 
set U of Aid, we have 


lim lim 

T —>-oo N—*oo 


1 

TN^ 




IT 


-W: 


N 


G C 


< 


lim lim 
T —^00 N —*oo 


1 

tW 


logP^w 


IT ^ 


G 


u 


— inf $(J| 7 ) , 
jgC 

> - inf 4 >(J| 7 ) . 
Jeu 








16 


L. BERTINI, A. DE SOLE, D. GABRIELLI, G. JONA-LASINIO, C. LANDIM 


If T were fixed, the above large deviation principle would directly follow from 
Theorem o The asymptotic T —> oo is related to the so-called T-convergence 
of the rate functions. We first discuss this issue in a general setting. Let T" be a 
metric space, we recall that a sequence Ft '■ X ^ [0, -l-oo] of functions T-converges 
to F ■. X [0, -l-oo] as r ^ oo iff for each x € X the following holds 

for any sequence Xt ^ x we have F{x) < lim Ft(xt) (4-3) 

T — 

there exists a sequence xt ^ x such that F{x) > lim Ft{xt) (4.4) 

T —>-oo 

Lemma 4.3. Let Pn,t be a two parameter family of probabilities on X endowed 
with its Borel a-algebra. Assume that for each fixed T > 0 the family {TW.tIatgn 
satisfies the weak large deviation principle with rate function TFt, that is for each 
K compact and U open in X we have 

lim ^logPM t{K) < —T inf Ft{x) (4.5) 

N^OO N ’ xGK 

lim log Pat t(C^) > —T inf Ft{x) (4.6) 

N^oo X ' xdU 

Assume also that the sequence Ft T-converges to F as T —>■ oo. Then for each K 
compact and U open in X we have 

lim lim ^ logPAr,T(i^) < - inf P(a:) (4.7) 

T— ^ooN^oc) I\ 1 x^K 

lim lim ^logPAr.T(C/) > - inf P(a;) (4.8) 

T^oo AT^oo XI xGU 

Proof. To deduce (lOi- pini from we need to show that for each K 

compact and U open in X we have 

lim inf Ft{x) > inf F{x) 

x^K 

lim inf Ft(x) < inf F(x) 

T^ocxGU xeu 

These bounds are a direct consequence of (IO)-(IOIl , see e.g. [Sj Prop. 1.18]. □ 


The following Lemma follows from Theorem by contraction principle. 


Lemma 4.4. Consider a profile 7 G C'^(T‘*) bounded away from 0 and 1 and a 
sequence {rj^ ■ N > 1} such that {rj^) j in F. Then for each T > 0 we have 
the following large deviation principle for the mean empirical current 'W^/T. For 
each C closed and each U open in A4d 


lim —7 log P: 




1 T 


G C 


< —T inf ^^(Jly) 
jgc 


lA 


A 


V — 




1 T 


G U 


> —T inf $t(J| 7) 
Je(7 


where we recall the functional ^t is defined in dnj. 


Proposition 4.5. Let 7 G C‘^{T‘^) bounded away from 0 and 1. The sequence of 
functionals 4)t(-|7) T-converges to the functional $(-|7) defined after 114.211 . 
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The previous Proposition, together with T;eTnmata J4.,Sl a,nrl 14.41 proves the “large 
deviations principle” stated in Theorem oi for compact sets. For its proof, we 
need a few preliminary Lemmata. For each 3 ^ B, since the empirical density is 
bounded, by the continuity equation 12.411 . we have that $t(J| 7 ) = +oo if T is 
sufficiently large. We next show that this holds uniformly for all J whose distance 
from B is uniformly bounded below. To this end we introduce the following metric 
on Md- Pick a sequence of smooth vector fields Gk & (7^ (T'^; , fc > 1, dense in 

the unit ball of for J, J' S Aid we then define 

CX3 ^ 

0 (J,J') = ^^lA|(J-J',Gfe)| 

fe=i ^ 

It is easy to show that p is a metric inducing the weak topology of Aid- 

Lemma 4.6. For each 6 £ (0,1) there exists Tq = Tq(S) £ ]R._|_ such that for any 
T > Tq we have <i>T(J| 7 ) = +oo for any 3 £ Aid such that g{3,B) > 6. 

Proof. Given J £ Aid let us denote by J £ S its projection on the subspace B. If 
3{du) = jdu for some j £ this is simply the orthogonal projection of 

j to B. In general J £ Aid is defined by (J,VF) = 0 for any F £ (7^(T‘^) and 
(J,E) = (J,E) for any E £ such that V • E = 0. It is easy to verify 

that J is uniquely defined by the above requirements. We then have 

oo .. oo 

(5 < g{J,B) < e(J, J) = E ^ 1 ^ |(J - J,Gfc)| = ^ ^ 1 A |(J, VFfc)| 

k=l k=l 

where F^ £ (T*^) is obtained from Gk by solving the Poisson equation AFk = 

V-Gfc so that Gk = VFfe+Efc with V-Ej, = 0. Note that \\Fk\\L 2 < Go||G/c||l 2 < Cq 
for some constant Cq not depending on k. 

From the previous inequality we get that there exists k = fc(J) such that 
|(J,VF^)| > 5. Let W £ 71t,j and denote by tt*, t £ [0,T] the corresponding 
solution of the continuity equation itnil . By choosing F = Fj. and using that 
Wt=T3, from (EJ) we get 

{7TT-^,Fk)=T{3,VFk) 

Since —I < ttt — J A !> the absolute value of the l.h.s. above is bounded above by 
Cq. On the other hand, since g{3,B) > 6, the absolute value of the r.h.s. above is 
bounded below by 6T. By taking Tq > Cq S~^, the Lemma follows. □ 

Lemma 4.7. Consider a profile 7 £ C'^(T‘*) bounded away from 0 and 1 and 
let 3 £ B, T > 0. We then have $t(J| 7 ) < +00 iff 3{du) = jdu for some 
j £ L 2 (T‘^;R‘^). Moreover there exists a constant Ci £ (0, 00 ) (depending on 
such that for any T > 0 and any 3{du) = jdu we have 

7^(jJ)<<i>T(J|7)<t^i[(jJ) + l] (4.9) 

£1 

Proof. Let W £ At.j', by choosing the vector held F in the variation expression 
EH) constant in time and divergence free we get 

i/[o,T](W|7)>(J,F)-ii(F,F) 
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where we used that x( 7 r) <1/4. Recalling that J S S, by optimizing over F S S we 
see that $t(J| 7 ) = +cx) unless 3{du) = j c?u for some j S L 2 (T‘^;K‘^). In fact this 
argument also proves the first inequality in 14.911 . To prove the second inequality 
in (lOli it is enough to construct an appropriate path W G At, 3 ', we simply take 
Wtidu) = tjdu. The solution of the continuity equation (12.411 is then given by 
TTt = 7 and, by Hd. 6 |l . 


i/[o,T](W|7) = i( ^ 


X(7)’ 


j + jV, 


Recalling that 7 is bounded away from 0 and 1, xil) = 7(1 “ 7)j the result follows. 

□ 


We next show that on divergence free measures J the functional T$ 7 ’(-| 7 ) is 
sub-additive. 


Lemma 4.8. Consider a profile 7 G C'^(T‘^) bounded away from 0 and 1 and let 
3 G B. Then, for each T,S > 0, we have 

(T + .5) $t+s(J| 7) < T$r(J|7) + ^$s(J|7) 

Proof. By (14.111 . given e > 0 there are G At ,3 and such that 

4>t(J| 7) > r/[o.T](W|7) - ‘&s(J|7) > S 

Let Wt, t G [0,T -I- S'] the path obtained by gluing with W^, i.e. we define 
Wt := + lI[T,T-i-S](^)W/_jn, and denote by the corresponding solution of 

the continuity equation dm . Then W G .4.t+s,j and, by the invariance of d[o,T] 
with respect to time shifts, 

(r + S)$r+s(J|7) < /[o.r+s](W|7) = J[o,t](W1|7)+/[o,s](WVt) 

< r$T(J|7) + ^$s(J|7)+e 

where we used that 3 G B, G At, 3 implies ttt = 7 . □ 

Recall that 4 )(-| 7 ) is defined in 14.211 and note that, by Lemmata 14.til and K7\ 
<i)(J| 7 ) equals -l-oo unless 3{du) = jdu for some divergence free j G L 2 {T‘^; By 
the sub-additivity proven above we have that <i>(J| 7 ) = limT—joo 4 ’t(J| 7) pointwise 
in J. However the pointwise convergences does not imply the T-convergence and 
some more efforts are required. 


Lemma 4.9. Consider a profile 7 G C^(T‘^) bounded away from 0 and 1, and let 
3 G B. Then, for each open neighborhood U of 3 we have 

lim inf $t(J| 7 ) < inf <i>(J| 7 ) 

T^oo Je(7 3eu 


Proof. Thanks to Lemma ^21 we can assume that the set 17 is a bounded subset of 
^ 2 ( 1 ’'^;®'^)- Pick S > 0 and let k := [T/S], R := T — kS G [0,S); by Lemma f4.81 
we get 


lim inf <i)T(J| 7 ) < 
T^oo Jg£/ 


lim inf 

fe^ooJe(7ne 1A:S-|-R 


kS 


$s(J|7) 


R 

kS + R 


‘i’R(J|7) 


< 


lim 

k —>-oo 


kS 


inf 


IfcS -I- i? Jec/nB 


^s(J|7) + 


R 

kS + R 


sup 4 'r(J| 7 ) 
JGC/nB 


= inf 4 * 5 ( 417 ) 
3eunB 
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where we used again Lemma W7\ to get supjgjy^g <i>_R(J| 7 ) < oo. By taking the 
infimum over S' > 0 we get the result. □ 

Lemma 4.10. Given m € (0,1), let Sm '■ K_|_ be the functional 

Sm{p) ■■= f du |p(u)log^^ + [1 - p(u)] logi— 

Then for each i5 > 0 there exists Tq = Tq(6) > 0 such that the following holds. For 
each 7 i ,72 G iFm there exists a path (W, tt) £ C([Q,Tf\\ M.d x T) n2l.y^, such that: 
TTO = 7l> ^To = 72, 

|(Wt„,F}| < iro||V-F|U,+5||F|U, foranyF€C\T^-X), (4.10) 

and 

VTo]((W,7r)|7i) <S™(72)+5 (4.11) 

Proof. The strategy to construct the path tt is the following. Starting from 71 we 
follow the hydrodynamic equation 12 . dt until we reach a small neighborhood (in a 
strong topology) of the constant profile m, paying no cost, then we move “straight”, 
paying only a small cost, to a suitable point in that small neighborhood which is 
chosen so that starting from it we can follow the time reversed hydrodynamic 
equation to get to 72 ; the cost of this portion of the path is Sm(,l 2 )- The current 
W is chosen so (im . hold, i.e. is the one whose cost is minimal among the 
ones compatible with the density path tt. 

Let A £ Tm and denote by PtX the solution of the Cauchy problem 12.311 {Pt 
is indeed the heat semigroup on T"^). By the regularizing properties of the heat 
semigroup, given > 0 there exists a time Ti such that ||PtA||ff^ + ||PtA — m||oo < 
(5i, for any t > Ti. Here = ||Vi^||l 2 is the standard Sobolev norm on 

T''* and the time Ti is independent on A because 0 < A < 1. We now choose 
(5i < (1/2) [m A (1 — m)] and let To := 2ri + 1, 7 ^ = PtiH, j = 1, 2. The density 
path TT is then constructed as 

[ Pt7i for t £ [0, Pi] 

TTt := i 7 i[l - (t - Pi)] + 72 [t - Pi] for t£(Pi,Pi + l) 

[PTo-t72 for t£[Po-Pi,Po] 

while the associated current path W is such that 

for t£[0,Pi] 

Wt= l-XVTTt+x{7 Tt)^Ht for t£(Pi,Pi + l) 

[ iVPTo-t 72 for t £ [Po - Pi, Po] 

where H £ + 1) x T'^) solves 

V • [x(7rt)VPt] = -dt^t + ^Att* 

Note there exists a unique solution since the r.h.s. is orthogonal to the constants 
(note 7 i £ Frn, i = 1 , 2 ). 

It is straightforward to verify that (tt, W) £ 21.^^ , i.e. the continuity equation 
(EH holds. Thanks to the invariance of I with respect to time shifts, the cost of 
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this path is 

ho,To] ((W,^)| 7i) = -f[0,Ti]((W,7r)|7i)+/[0,i]((W._Ti,7r._ri)|7i) 

+ ^[o,Ti]((W._(Ti+i),7r._(rj+i))|72) (4-12) 

Since in the time interval [0, Ti] the path follows the hydrodynamic equation, the 
first term on the r.h.s. of vanishes. By considering the time reversal of the 

portion of the path in the interval [Tq — Ti, Tq], it is straightforward to verify, see 
Pl Lemma 5.4], that 

4[o.ri]((W._(rj+i),7r._(T^_,.i))|72) = £'^( 72 ) - £m(72) < £m(72) 

It remains to bound the second term on the r.h.s. of 14.1211 . Note that, by 
construction, in the interval (Ti,Ti + 1 ) we have 


inf 7 i(m) a inf 72 (u) < ttj < sup 71 (m) V sup 72 (u) 

U U q, ,, 


which implies, by the choice of i5i, 

i [m A (1 — to)] < TTt < 1 — i [to a (1 — to)] 


i.e. the the path is uniformly bounded away from 0 and 1. By the same compu¬ 
tations as in Pl Lemma 5.7], it is not difficult to show that there exists a constant 
C > 0 depending only on to such that 


4[o.i] (W._Ti,7r._Ti|7i) 


/■i X 

< C dt \\dtTTt+Ti + 2^^t+Ti 

< c'[||7.||f+IAA] 


||2 


which, by taking (5i small enough, concludes the proof of lUTTll . The bound (imH) 
follows easily from the construction of the path W by using Cauchy-Schwartz and 
what proven above. □ 


Lemma 4.11. Consider a profile 7 G C^(T^) bounded away from 0 and 1, and let 
3 G B. Then, for each open neighborhood U of 3 we have 

lim inf 4 >t(J| 7) > inf 4 >(J| 7 ) 

JGC/ 

Proof. Recalling definition (lOl . it is enough to show 

lim inf $t(J| 7 ) > lim inf $j’(J| 7 ) (4.13) 

T^oo-I6(7 

Given T > 0 there exist 3^ G U and G At.j^ such that 
m^4>T(J|7)>^/[o,r](Wi|7)-l 

Let TTj, t G [0, T] the density path associated to W^, 5 > 0 and To as in Lemma. li.l 01 
We now define, on the time interval [0, T -I- To], the path 

Wt := + I[T.T+T„](i)W?_T 

where W^, t G [0,ro] is the path constructed in Lemma OUI with 7 i = TT^ and 
72 = 7- Let finally J = Wt+To/(T + To) = TJV(T + Tq) + W%/iT + Tq); note 
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that, by construction, J G B. From T;cmma, f4.1 01 it now follows that, for T large 
enough, 3 G U and 

-^/[o,T+To](W|7) < + -^-^[SM + S] ■ 

1 lo i + io i + io 

By taking the limit T ^ oo, follows. □ 


Proofs of Proposition \J^ and Remark \4-1\ By Lemmata 14. til and 14. Ill we have 

that, for each J G Aid and any neighborhood 17 of J 

lim inf $t(J| 7 ) = inf $(J| 7 ) 

T^oo J6(7 JeU 


The F-convergence of the sequence $t('| 7 ) to the lower semi-continuous envelope 
of d)(-| 7 ) now follows from the topological definition of F-convergence, see e.g. |S1 

§1-4]. 

We next prove Remark IO Let m G (0,1); by using the path introduced in 
Lemma [4. 101 it is straightforward to check that, for each 71,72 G Am and J S we 
have $(J| 7 i) = $(J| 72 ) which proves the first statement. 

Since d>(-| 7 ) is the lower semi-continuous envelope of <I>(-| 7 ), it is enough to prove 
the convexity of the latter. As is a closed convex subset of Aid, it is furthermore 
enough to show that for each Ji, J 2 G B, each p G (0,1), and each 7 £ C'^(T‘^) 
bounded away from 0 and 1 , we have 

$(p Ji + (1 -p) J 2 I 7 ) < p 5 (Ji| 7 ) -I- (1 -p) 5 (J 2 | 7 ) (4-14) 

Given e > 0 we can find T > 0, G ApT.Ji, and £ -^{i-p)T ,32 so that 

5(Ji| 7) > ^4[0,pT](Wi|7)-e 

5(J2|7) > (1 4[o,(i-p)T](W^l7) - £ 

By the same arguments used in Lemma [4.81 the path obtained by gluing with 
is in the set 7lT,pji-i-(i-p) Ja- The bound (14.1411 follows. □ 

We conclude this section proving the exponential tightness needed to complete 
the proof of Theorem 10 


Lemma 4.12. Fix a sequence 
{Ki : ^ > 1 } of Aid such that 

lim lim - 7 

T^oo N^oo TN’^ 


G Xn- 


logP^f^ 


There exists a sequence of compact sets 




< -£■ 


Proof. Fix a vector field H : ^ We claim that for every A > 0, 


P: 


|(W. 


^,H)| > AT 


< 2exp{-TN‘^[AC{U)-^-7]} 


(4.15) 


where (7(11) = max^jl V ||77j||^ V \\dujHj\\^}. Indeed, by Chebychev exponential 
inequality, 


P: 


,N 


(W^,H) > AT 


^ -SATN'^TffN 

S e 


A4T(6',H)e^"(Tp,H) 
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for every 6* > 0. Here, is the mean one exponential martingale defined 

just after (ESI), with a now time independent vector field H, and 

d 


f=l xST^ ° 


t=i xeTf, 


A Taylor expansion shows that the absolute value of / 7 v(T, r^jH) is bounded by 
C(li)TN‘^{9 + 29^e^^^^^}. To conclude the proof of the claim it remains to choose 

9 = to remind that Mt{9,'H.) has mean one and to repeat the same 

argument with — H in place of H. 

Recall the definition of the sequence {Gfc : fc > 1} defined just after Proposition 

10 and assume, without loss of generality, that C'(Gfc) < k. For each £ > 1, the 
set Ki of measures defined by 

Ke = f|{J : |(J,Gfe)| < ik + 7f£} 

k>l 


is compact. On the other hand, by 14.15II . 

1 


t,N 
■ r,N 


IT 


W;^ G ATf 


< 4e 


-TN'^i 


provided N is sufficiently large. This proves the Lemma. 


□ 


5. Dynamical phase transitions 

In this Section we analyze the variational problem defining the functional 
$. For the symmetric simple exclusion process we prove, in Subsection o that 
$ = U, where U is defined in EH . Therefore no dynamical phase transition 
occurs in this case. In Subsection E3 we consider a system with general transport 
coefficients and show that, under suitable convexity assumptions, it is possible to 
construct a traveling wave whose cost is, for J large, strictly less than the constant 
profile. These convexity hypotheses are satisfied for the KMP model EEl and 
therefore we prove that a dynamical phase transition takes place. 

5.1. Symmetric simple exclusion process. The following statement is essen¬ 
tially proven in for the reader’s convenience we reproduce below its proof in 
a more formal setting. Together with Theorem 14.21 concludes the proof of Theo¬ 
rem |0| 

Proposition 5.1. For each m G (0,1) the functional Um '■ M.d —> [0, cxd] defined 
in (123 is lower semi-continuous. Moreover if j H Tm is bounded away 

from 0 and 1 we have Um{3) = ‘h(J|7) for any J G M-d- 

Proof. We first prove the lower semi-continuity of Um- Given (J,p) G Md x Pm 
we define 

:= supVf(p,J) 

F 

Vf(p,J) := (J,F)-i(p,V-F)-i(F,x(p)F) 


where 
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and the supremum is carried over all smooth vector fields F G Note 

that, if 3{du) = jdu for some j G L 2 (T‘^;R'^) and p G C^{T^) is bounded away from 
0 and 1 , we have 

«(/-■ J) = i ([j + t Vp], [J + i Vp]) 

Recalling definition (ira . by the approximation arguments in Subsection roi we 
have t7r„(J) = infpg^^ J). 

By the concavity of x{p) we have that, for each fixed F, the functional Vf(-, •) : 
T X M-d —> R is convex and lower semi-continuous. The lower semi-continuity of 
hence of Um, follows now easily. We also note that the previous argument shows 
that U is a convex functional on !F x Aid- 

We next prove that for each 7 G C'^(T'^) D J-m bounded away from 0 and 1 we 
have f7m(J) = ‘h(J|7)- From the definitions 12.811 . 14.211 and Lemma \n\ we can 
assume that 3{du) =idu for some j G L 2 (T‘*;R^) divergence free. 

We first show that for each T > 0, and each path (W, tt) such that W G At.j 
we have 

i/[o.T]((W,7r)|7) > C/^(J) 

Indeed, thanks to the approximation constructed in Section fit. 41 we can assume that 
TT is a smooth path bounded away from 0 and 1. For such a smooth path (ESI) yields 

i/[o.T] ((W, T)|7) = 1 £dt Ui-Kt, Wt) >U(^^ £dt Um{J) 

where we used, in the second step, the joint convexity of the functional U and 
Jensen inequality. In the last step we finally used that, by conservation of mass 
TTt G Am for any t G [0, T], 

To show the converse inequality it is enough to construct, for each T large enough, 
an appropriate path. Given e > 0 there exists p G C'^(T‘^) D Am bounded away 
from 0 and 1 such that Um{jdu) > ZJ(p,j) — e. For T > 2 we construct the path 
(W, tt) such that 

{ •wdu if t G [0,1) 

^idu if tG[l,r-l] 

—wdu if t G {T — 1,T] 

where w solves V • w = 7 — p. It exists because 7 , p G Am, i.e. they have the same 
mass. The density path tt is the corresponding solution of fIM . i.e. 

( 7(1 -t)+pt if t G [0,1) 

TTt = i p if t G [1,T - 1] 

[ p{T-t)+x{T+l-t) if tG{T-l,T] 

It is straightforward to verify that W G At,j- Moreover 

hm i/[o_^j((W,^)| 7 ) =U{p,3) < Um{J)+e 

1 —^00 J. 

which concludes the proof. □ 

We conclude this Section by showing that in the one dimensional case the func¬ 
tional Um is given by EiS 
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Lemma 5.2. Let m G (0, 1), d = \ and J{du) = j du for some j G K. Then 

1 f 

Um{J) = 


2xim) 


Remark. As it will be apparent from the proof, this Lemma holds whenever the 
real function p i—> l/x{p) is convex. 


Proof. Let p G C^(T) n be bounded away from 0 and 1. We have 


du 


[j + {l/2)p\u)]^ 

X{p{u)) 


= du , . + du 

Jt xkpw) Jt 


X{p{u)) 


because the cross term vanishes upon integration. By Jensen inequality, 


X{p{u)) 


> 


J 


X(m) 


On the other hand, by considering the constant profile p{u) = to, we trivially have 
UmU du) < (1/2) j^/x(TO). The lemma is therefore proven. □ 


5.2. Other models. The general structure of the hydrodynamic equation obtained 
for the scaling limit of the empirical density for stochastic lattice gases with a weak 
external field E has the form, see uniiis] 

dtp + V- W{p) = 0 

W(p) = -izJ(p)Vp + x(p)E 

where D(p) is the diffusion coefficient and xip) is the mobility. 

In this general context, for smooth profiles, we let 

U{p, J) = i ([J - W(p)], ^ [J - W(p)] ) 

The integrated empirical current is expected to satisfy, see ^ for a heuristic deriva¬ 
tion, a large deviation principle with rate function 

Vt](W)= [ dtU{TTt,Wt) 

Jo 

in which tt is obtained from Wj by solving the continuity equation dtTT + V • W = 0. 


We analyze the variational problem (lOl in this general setting and show that, 
under some assumptions on D{p),xip), a time dependent strategy is more conve¬ 
nient than taking a density path tt constant in time, so that <!><[/. For simplicity, 
we here discuss only the one dimensional case d = 1 and assume that there is no 
external field, E = 0; see El § 6.2] for more details. 

Given a mass to and u G K., let ^ be defined by 

{J + v [p{u) - to ] - w{p{u))y 


^viJ) = 


inf i 

p 2 . 


du ■ 


X{p{u)) 


(5.1) 


where w{p) = W{p) = —{l/2)D{p)yp + x{p)E and the infimum is carried over the 
profiles p of mass to, i.e. over Fm- It will be convenient to write the term D{p)\Ip 
as Vd{p), i.e. d{p) is an antiderivative of D{p). 

We claim that for each u G M 


$ < J/ 


(5.2) 
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Indeed, consider a profile po in Tm- Let T = v~^ and set p{t,u) = po(u — vt), 
w(t, u) = J+v[po{u—tv)—m\ in the time interval [0, T], An elementary computation 
shows that the continuity equation holds and that the time average over the time 
interval [0,T] of w{-,u) is equal J. In particular, 

^(J) < ^ dtU{p{t),w{t)) . 

On the other hand, it is easy to show by periodicity that the right hand side is 
equal to 

1 {J + u[po(u) - m] - w{po)y 

2^0 ^ x(Po(m)) 

Optimizing over the profile po, we conclude the proof of itCTl . 


We next show that, if the real function p l/x{p) is convex and }("{rn) > 0 for 
some 0 < TO < 1, then $(J) < U{J) for \ J\ sufficiently large. 

To prove the previous statement, we first note that, in view of 115.211 and of 
Lemma 15.21 it is enough to show that there exists A £ M such that 


lim sup 

I J|—»'C!0 


J2 


1 

2x(to) 


(5.3) 


Fix a mass to, a current J and take v = XJ. For p £ Trm by expanding the 
square we get that 


{■/ + AJ[ p - to] + (l/2)Vd(p)}2 


= r 


du 


x(p) 

{1 + A[p-to]}2 1 

X{p) 4 


(5.4) 


du 


mp)f 

x{p) 


because the cross term vanishes. Expand the square on the first integral. Let 
F{r) = be the smooth function defined by 


^ {l + A[r-TO]p _ 

^ ^ X{r) 

An elementary computation shows that 

F”im) = --^^|2x(to)^A^ - 4x(to)x'(to)A + 2x'(to)^ - x(w)x"(w)| . 


Let A = X^(™)/x(^)- For this choice F"{m) < 0. In particular, we can choose 
a non constant profile p{u) in Tm close to to such that F"{p{u)) < 0 for every 
u. Hence, by Jensen inequality, the coefficient of in 15.411 is strictly less than 
x(to)“^. This completes the proof of the claim. 


For the KMP model 0^3 have D{p) = 1 and x(p) = P^■ In particular x 
and 1/x are convex functions. Hence, by the above results, we have 4)(J) < U{J) 
for all sufficiently large currents J. 

For the weakly asymmetric exclusion process with large external field a similar 
phemenon occurs. More precisely, as shown in |7], there exists a traveling wave 
whose cost is strictly less than the constant (in space and time) profile. A numerical 
computation [Jj suggests also that the minimizer of the variational problem (lOl 
is indeed a traveling wave. 
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